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Q ; ABSTRACT 

o 
o 

Beamed 7-ray burst (GRB) sources produce a bow shock in their gaseous 
Q\ ', environment. The emitted flux from this bow shock may dominate over the direct 

C^ I emission from the jet for hues of sight which are outside the angular radius of the 

0> I jet emission, 0^. The event rate for these hues of sight is increased by a factor of 

~ 260(Ob/5°)~'^. For typical GRB parameters, we find that the bow shock emission 
from a jet of half-angle ~ 5° is visible out to tens of Mpc in the radio and hundreds of 
Mpc in the X-rays. If GRBs are linked to supernovae, studies of pecuhar supernovae 
^D . in the local universe should reveal this non-thermal bow shock emission for weeks to 

months following the explosion. 
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^ : 1. INTRODUCTION 



The afterglows of Gamma-ray bursts (GRBs) are most naturally described by the so-called 
"fireball" model (see e.g., Paczyhski & Rhoads 1993; Katz 1994; Meszaros & Rees 1993, 1997; 
Waxman 1997a, b; Sari, Piran, & Narayan 1998). In this model, a compact source releases a large 
5^ I amount of energy over a short time and produces a relativistically expanding fireball. Eventually 

the fireball interacts with the circumburst medium, producing a spherical relativistic shock in 
it. As the shock decelerates due to the accumulation of mass from the external medium, it 
approaches a self-similar solution (Blandford & McKee 1976) and produces delayed synchrotron 
emission in X-rays, optical and radio, similar to the observed afterglows. The model generically 
yields power-law spectra and lightcurves in general agreement with observations at X-ray (Costa 
et al. 1997), optical (van Paradijs et al. 1997) and radio (Frail et al. 1997) wavelengths. Precise 
positions have allowed redshifts to be measured for a number of GRBs (Metzger et al. 1997), 
providing a definitive proof of their cosmological origin. 

Recent observational evidence indicates that at least some GRB are not spherical explosions. 
A broad-band break in the lightcurve power-law index was predicted for shocks produced by 
collimated jets due to the lateral expansion of the jet when its Lorentz factor drops below the 
inverse of its opening angle (Rhoads 1997, 1999a, b; Panaitescu &: Meszaros 1999). Such breaks 



have been seen in GRB 990510 (Stanek et al. 1999; Harrison et al. 1999), GRB 991216 (Halpern 
et al. 2000), and GRB 000301C (Sagar et al. 2000; Masetti et al. 2000; Jensen et al. 2000; Berger 
et al. 2000). The ratio of the spectral index to the light curve index also suggests non-spherical 
energy ejection for some events (e.g. GRB 991216: Garnavich et al. 2000). 

The possibility that GRBs are associated with a rare sub-class of supernovae (SNe) was 
advocated recently based on a possible SN component in the lightcurves of GRB 980326 (Bloom 
et al. 1999) and GRB 970228 (Reichart 1999) and the association of SN1998bw with GRB 980425 
(Galama et al. 1999; Kulkarni et al. 2000). In this case the GRB emission needs to be beamed in 
order to achieve the high Lorentz factor required for GRB outflows, despite the massive progenitor 
envelope and the low energy yield of SNe (MacFadyen, Woosley, & Heger 1999; Khokhlov at al. 
1999; Umeda 2000). 

The current state of knowledge leaves two major questions open: (i) are GRBs associated 
with rare SNe?; and (ii) if so, what is the characteristic collimation angle of the energy release? 
The collimation angle has important implications on the nature of the central engine. If the 
average angular radius of the two opposing jets in collimated GRBs is Qb, then the total energy 
release will be reduced by a factor 2tt@1/4:Tt = 0^/2 and the event rate will be increased by the 
inverse of this factor. In this paper we propose a direct observational probe that will help answer 
these questions. 

When a transient relativistic jet (equivalent to a relativistic "bullet") moves into the 
circumburst medium (be it the interstellar medium or a progenitor wind), a bow shock is generated 
in front of the jet. The bow shock accelerates electrons which emit synchrotron radiation, similarly 
to the primary shock in front of the jet. For observers situated within the cone of the jet emission, 
the bow shock emission will be sub-dominant relative to the jet emission due to relativistic motion 
of the jet towards the observer. However, for highly collimated outflows, most lines of sight lie 
outside the jet cone, and for those - the bow shock emission may dominate. In particular, if GRBs 
originate in SNe then the bow shock emission will add a strong non-thermal component to the 
emission from the SN remnants extending from the radio to the X-ray bands. This emission will 
extend over much longer times compared to the jet emission since the Lorentz factor of the bow 
shock is typically much smaller than that of the jet, resulting in a smaller relativistic compression 
of time in the observer's frame. In the following sections we will apply the same model used 
ordinarily for the primary fireball shock, to calculate the flux in this non-thermal bow-shock 
component. In §2 we present our model for the bow shock emission; in §3 we show our numerical 
results; and finally, in §4, we summarize our main conclusions. 



2. MODEL 

The impulsive ejection of a relativistic jet from the GRB source results in a thin disk of 
material moving outwards and expanding, which we refer to as a relativistic bullet. As long as 



the angular radius of the bullet relative to the center of the explosion, Qb, is much larger than 
the inverse of its Lorentz factor, the bullet will expand as if it is part of a spherically symmetric 
fireball. This follows from the relativistic transformation of velocities, which implies that the 
propagation speed of a signal in the perpendicular direction to the bullet motion cannot exceed 
~ c/'ji) in the frame of the stationary ambient medium. For simplicity, we limit our discussion to 
this regime. 

In our model we assume that as the relativistic bullet moves through the ambient medium, 
it deposits a fraction of its kinetic energy in each infinitesimal distance it travels. We treat this 
energy deposition as a sequence of point explosions, which drive shocks into the surrounding 
medium. The combination of this sequence of small shocks results in a conical bow shock structure 
behind the bullet. 

Figure 1 illustrates the geometry of the bow shock. We use cylindrical coordinates, with the 
z-axis connecting the bullet center to the source. The distance between the source and the bullet 
is denoted as Rh- We assume that observers are located perpendicular to the 2:-axis, with 9 = 
along the source-observer axis. Next we consider a slice of the bow shock geometry perpendicular 
to the z-axis. The slice is a cylindrical shell of height 6z, thickness A and inner radius r. As the 
bullet passes through a point z it has a radius l{z) and a Lorentz factor 7b(z). The amount of 
thermal energy 5E which gets deposited by the bullet in the surrounding medium as it moves an 
extra distance 5z, is then given by 

SE = (75 — 1)-itI 6z rifiefnpC , (1) 

where n is the electron number density of the external medium, and ^e^T-p is the ion mass per 
electron. We assume that a fraction e of the energy stored in the "causally-connected" edge of 
the bullet is transferred to the external medium. The "causally-connected" edge of the bullet is 
defined to be the region occupying an angular size I/75 from the bullet edge. The energy deposited 
in the infinitesimal distance interval 6z drives a shock into the ambient gas. The gas behind the 
shock is piled into an outgoing thin shell. The shock expands perpendicular to the z-axis, but due 
to the time delay between the energy deposition at different points along the z-axis, the sum of the 
fronts generated by these points defines a conical shape for the resulting bow shock. We denote 
the Lorentz factor of the outgoing shell in the frame of the unshocked gas by 7 . The electron 
number density n' and internal energy density e' of the shocked gas in the frame comoving with 
the shell can be written as (Blandford & McKee 1976) 

, 77+1 fr,-, 

n = -n, (2 

7-1 

e' = — -(7 - l)nfiempc'^, (3) 

7-1 

where 7 is the adiabatic index of the shocked gas which changes from 4/3 and 5/3 as the shock 
velocity changes from the relativistic to the non-relativistic regime. We interpolate between these 



regimes using the simplified expression 7 « (47 + l)/(37) (Dai, Huang & Lu 1999). Equations (|2|) 
and (^ then become 

n' = i'jn, (4) 

e' = 467 rifiefnpC , (5) 

where 5 = (7 — l)/7. Note that equations (0) and (|5|) are appropriate for both the relativistic and 
non-relativistic regimes. The total kinetic energy of the shell is 

Ek = vr(r —l)5z rifie'mpC /3 7 , (6) 

where P"^ = 1 — 1/7^. Conservation of energy implies, Ek = e^6E, where 

29fc7fc - 1 . . 1 .„. 

is the ratio between the volumes of the causaly-connected region near the edge of the bullet (with 
angular size l/7b) and the entire volume of the bullet. We therefore obtain 



7 = ^1 + ^^7. - l)^, = y 1 + f:;^,, (8) 

where / = e^(7b — !)• The thickness of the shell is derived by using the conservation of particle 
number. This yields 

Figure 2 illustrates the geometry of the emission from an infinitesimal volume element in 
cylindrical coordinates, dV = rdrdOdz. We define the emission coefficient j'^, to be the power 
emitted per unit frequency, z/', per unit volume per steradian in the rest frame of the emitting 
material. We use prime to denote quantities in the local rest frame of the emitting material, while 
unprimed quantities are measured in the rest frame of the external medium. Note that jy/v"^ is 
Lorentz invariant (Rybicki & Lightman 1979). For each slice, the cylindrical shell of expanding 
material emits isotropically in its local rest frame with f^, = P'(z^', r,t) /47r and v' = 1^7(1 — /3/u), 
where 7 and f3c are the Lorentz factor and the velocity of the emitting material, and /i = cos 9. A 
photon emitted at time t and place r in the unshocked gas frame will reach the detector at a time 
T given by 

T. = :r^ = *-^. (10) 

1 + z c 

where Z is the cosmological redshift of the GRB and T is chosen so that a photon emitted at the 

GRB source at t = will arrive to the detector at T = 0. Thus we have |^ 



Note that a similar equation for calculating the emission from GRB afterglows in a spherical coordinate system 
was originally derived by Granot, Piran, & Sari (1999). 



where D is the luminosity distance to the GRB, and 7, /3, /x are evaluated at the time t implied by 
equation (10). The factor of 2 is added to describe the combined emission from the bow shocks 



generated by two opposing jets. 

The volume integration expressed in equation (^) should be taken over the region occupied 
by the emitting bow shock at a given observed time. The lower and upper boundaries of the 
integration over r are derived by considering the relativistic time delay (similar to the approach 
used by Wang &; Loeb, 2000). For a given observed time T, the outer boundary Ri satisfies the 
equation 

T = -F{R^)-^, (12) 



c 



where 



The inner boundary R2 can be obtained by solving the equation 

Next we derive the local synchrotron emissivity in the two regimes where the infinitesimal 
segment of the bow shock under consideration is moving at relativistic or nonrelativistic speeds. 
In the relativistic regime, we assume that the energy densities of the shock-accelerated electrons 
and the magnetic field are fixed fractions of the internal energy density in the shell behind the 
bow shock, Cg = Eec', e'^ = e^e', and that the bow shock produces a power law distribution of 
accelerated electrons with a number density per Lorentz factor of A^(7e) = K^~^ for 7e > 7mini 
where 

7min=(^);^, K = {p-l)n'^^^l (15) 

When the bow shock is nonrelativistic, we use eg to denote the fraction of the internal energy 
which goes to relativistic electrons with 7 ;^ 2. Again, these relativistic electrons obtain a power 
law distribution A^rei(7e) = K'y~P for 7e > 7miin where 

rnp 



K = ee^/x,(p - 2)47(7 - l)n7^„', (16) 



and 7min = 2. Hence the number density of relativistic electrons is 

n ,1 = ee^^ugi-— 47(7 - 1) . (17) 

me p-l 7min 

This prescription describes well the accelerated electron population in the non-relativistic shocks 
of supernovae (Chevalier 1999). The characteristic values inferred for the parameters p ~ 2-3 
and Ce ~ 1-10% in supernova shocks (Chevalier 1999; Koyama et al. 1995, 1997; Tanimori et al. 
1998; Muraishi et al. 2000) are in the same range as those for GRB afterglows (Waxman 1997a, b; 



Wijers, Rees, & Meszaros 1997; Sari, Piran, &; Narayan al. 1998). We will consider different 
parameter values, but for simplicity, we will adopt the same values of eg and p in describing 
both the relativistic and sub-relativistic regimes. This is justified by the similarity between the 
parameter values which are typically chosen to describe relativistic GRB afterglows (e.g., Waxman 
1997a,b) and non-relativistic radio supernovae (e.g.. Chevalier 1998). 

At sufficiently high Lorentz factors, the electron cooling time is shorter than the dynamical 
time. The critical Lorentz factor 7c above which electrons cool on a time shorter than tdyn is given 
by (Sari, Piran & Narayan 1998) 

3me 1 , . 

IGeBaTfJ-empC UynOJ n 

where idyn is equal to the dynamical time in the frame of the external medium, i.e., the time it 
takes for a cylindrical shell to expand from an initial radius / to r, 

tdyn = ^F{r). (19) 

Note that equation (|l8|) is different from equation (6) in Sari et al. (1998) because of the 
introduction of b and /ig- An electron with an initial Lorentz factor 7e > 7c cools down to 7c in 
the time tdyn- 

The radiation power and the characteristic synchrotron frequency of a randomly oriented 
electron with a Lorentz factor 7e are given by (Rybicki &; Lightman 1979) 

P'ile) = \aTllP'c^, (20) 



IGrripC 



^syAle) = ,^_ ^ , (21) 



where B' = JSne'^, and me,qe are the electron mass and charge respectively. We define 

^c = ^syn(7c) ai^d ^'m = ^syn(7min)- In the regime of fast cooling, 'jmm > 7c, the emissivity is given 

by (Sari et al. 1998) 

K' = { Pl>,^.A^' /<)-'" K<^'<^'m (22) 






where 






(23) 



In the above equation, n' should be replaced by nj,j,j for the nonrelativistic regime. 
For slow cooling, 7c > 7min) the emissivity is given by 



Pi' = <! ^'.max W ly'mY'^'-^^''' ^'m < ^' < ^'c (24) 



At very low frequencies, synchrotron self-absorption becomes important. In the comoving 
frame of the shocked gas, the absorption coefficient ol^, scales as oi^, oc zy'^(P+^)/2 for v' > v'^ and 
as a^, oc v'~^l'^ for v' < u'^ (Waxman 1997a). We therefore write 



a',, = F^'-(^+4)/2, u'>u'^, (25) 

where (Rybicki & Lightman 1979) 

^^VM. f^S^Y^' im^cy-^KXB'i^^-^y^T ('J^) r ('J^±^) (26) 

Svrme V2vrm3c5y ^ " ^ V 12 / V 12 / ^ ^ 

and where A = (1/2) /g'^(sinQ)(*'+^)'^ sinada, and T{y) is the Gamma function. For u' < u'^ we 
then use 

al, = Hu'J^+'y^^) , u'<u'^. (27) 

Because va^ is Lorentz invariant, the absorption coefficient in the rest frame of the unshocked gas 
is 

a^ = 7(l-/3/i)a[,,. (28) 

Equation ( p!l| ) should then be modified to 

where Tj^ ~ a^ x A is the optical depth per unit frequency for synchrotron self-absorption across 
the shell thickness. 

The external medium could be either the interstellar medium (ISM) of the host galaxy 
(Waxman 1997a, b) or a precurser wind that was ejected by the GRB progenitor (Chevalier &; Li 
1999; 2000). For a wind profile p = Ar~'^, the distance between the source and the bullet is given 
by (Chevalier & Li 2000) 

Rb = 1.1 X 10^^ f — J E^2A~^ cm, (30) 

where 7?, is the bullet Lorentz factor, £'52 is the equivalent isotropic energy release of the GRB in 
units of 10^^ ergs, A = M„/AttV„ = 5 x 10^^^* g cm~^, Mw is the progenitor mass loss rate, and 
Vw is the wind velocity. A^, = 1 corresponds to M^ = 10~^ Mq yr~^ and V„ = 1000 km s~^. The 
electron number density is given by, n = p/{fiemp). Also note that the energy input into the two 
jets is equal to E x (G^/2). 

If the external medium is the ISM, Rf, is given by (Granot, Piran, &: Sari 1999) 

Rb = 5.53 X 10^^ (M^) '^' (^) '^' cm, (31) 

where ni is the ISM electron number density in units of 1 cm^^. 



3. NUMERICAL RESULTS 

We have solved the equations of §2 for different values of the free parameters in our model. 
In the following, we show results for a standard case, and five variations on the values of its 
parameters. In the standard case 1, we adopt the parameter values, @h = 0.1, £"52 = 1, e = 0.1, 
€e = 0.1, eB = 0.1. We also adopt A* = 1 and fi^ = 2 for the wind of a typical Wolf-Rayet star 
(Chevalier &: Li 2000). Cases 2-5 also refer to a wind profile, but with a variation on one parameter 
in each case compared to the standard case. In case 2 we adopt a lower density wind with A^, = 0.4 
(corresponding to V„ = 2500 km s~^ for the same value of M^ = 1 x 10~^ Mq yr~^). In case 3 we 
adopt, -E'52 = 10; in case 4, e = 0.5; and in case 5, p = 2.2. Finally, case 6 considers the ISM as the 
external medium with ni = 1 and yUg = Ij and all other parameters the same as in the standard 
case. 

In all cases, we start with 75 = 100, as appropriate for the afterglow phase. We show results 
only for the initial period during which 7^ ^ 1/06, since our model does not apply to later times 
when the lateral expansion of the bullet is important. Most of the initial kinetic energy of the 
bullet is dissipated in the ambient gas during this early period. In presenting our results we 
assume that the cosmological redshift of the source, Zg <C 1, because as it turns out, only nearby 
GRBs will produce sufficient bow-shock flux to be detectable. 

The numerical results for all six cases are shown in Figure 3. We plot the spectrum of the 
emission from the bow shock, with the vertical axis being the luminosity per unit frequency 
(Ljy = 4-kD^Fi^). The light curves appear to have similar features in all cases. At low frequencies, 
L,^ oc u"^ because of synchrotron self-absorption, while at high frequencies, L^, oc u^p'"^ due to 
efficient electron cooling. The primary peak in each lightcurve is located at the synchrotron 
self-absorption frequency (denoted as the peak frequency hereafter), while the second peak 
is located at the cooling frequency. Case 2 has a faster wind speed, hence a lower electron 
density than the standard case. Therefore, case 2 has lower peak luminosities and higher cooling 
frequencies. The observation time calculated in this case is also longer because it takes more 
time for the bullet to decelerate. Case 3 has a higher input energy than the standard case, but 
its peak luminosities are almost the same as those of the standard case. Also it has lower peak 
frequencies and higher cooling frequencies. This is caused by the fact that the bullet has more 
kinetic energy, and so its deceleration requires more time and a longer distance. As a result, 
the density of the wind is lower than that of the standard case when the bow shock emission is 
calculated. Case 4 assumes e = 0.5, implying that more energy is deposited into the bow shock 
by the bullet. Compared to the standard case, this case has a higher peak luminosity and a 
lower cooling frequency. Case 5 assumes a lower value for p. Compared to the standard case, it 
has a somewhat lower peak luminosity and a higher peak frequency. Case 6 considers the ISM 
as the ambient medium, implying a much lower electron density than the standard case with 
a wind profile. Consequently, the peak luminosity and the peak frequency are lower, while the 
cooling frequency is higher. Also note that in this case the peak luminosity increases with time, in 
contrast to the standard case. 



With these results at hand, we may now compare the thermal emission from the main 
supernova shock with the non-thermal emission from the GRB bow-shock in the optical regime. 
Type la supernovae (SNe la) have been used as "standard candles" because they all have similar 
light curves and peak absolute magnitudes. The absolute S-magnitude at maximum light (Mb) of 
typical SNe la is ~ —18.5 (Vaughan et al. 1995), corresponding to L^ ~ 1.3 x 10^® ergs s^^ Hz^^ 
at zv = 6.8 X 10^^ Hz. This is at least two orders of magnitude more luminous than the emission 
from the GRB bow-shock in our calculation. Type II supernovae can be generally divided into two 
relatively distinct sub-classes, SNe II-L ("linear") and SNe II-P ("plateau"). The majority of SNe 
II-L have a nearly uniform peak absolute magnitude, though there are a few exceptionally luminous 
SNe II-L (Young k Branch 1989; Gaskell 1992; Fihppenko 1997). The average Mb of the majority 
of SNe II-L is ~ —16.5 (Gaskell 1992), corresponding to L,y ~ 2 x 10^^ ergs s~^ Hz^^. This is still 
more than one order of magnitude brighter than the emission from the GRB bow-shock. But if the 
emission from the bow-shock persists for a sufficiently long time (as in case 3 of our calculation), 
then it may exceed the thermal SNe II-L emission due to the decline in the supernova lightcurve. 
SNe II-P show a very wide dispersion in the distribution of their peak absolute magnitudes; i.e., 
Mb spans the range from -14 to -20 (Young & Branch 1989), corresponding to L,^ ~ 2 x 10^^ 
- 5 X 10^^ ergs s^^ Hz^ . Thus, the less luminous SNe II-P (such as SN1987A), have peak 
luminosities in i?-band which are comparable to the emission from the GRB bow-shock. Of course, 
GRBs may occur in a rare subset of supernovae that have very different optical luminosities than 
the typical values mentioned above. 

The bow shock emission is more easily detectable at either radio or X-ray frequencies. At a 
frequency of ~ 5 GHz, most radio supernovae which reached their peak between 10 and 130 days 
had a peak spectral luminosity L^ between 10^^ and 10^^ ergs s~^ Hz~ (Li & Chevalier 1999, Fig. 
3). This peak luminosity is comparable to the bow shock emission in Figure 3 (except for case 5 
with p = 2.2, for which L,^ ~ 10^^ ergs s~^ Hz^ ). However, SN 1998bw, the most luminous radio 
supernova observed so far, reached after ~ 10 days a peak luminosity of L,^ ~ 10^^ ergs s^^ Hz^ 
at 5 GHz, which is about two order of magnitude brighter than our calculated bow-shock emission. 

In the X-ray regime, the best studied case of Type II supernova SN 1993J, had a luminosity 
of 3 X 10^^ ergs s^"*^ in the 0.1-2.4 KeV band on day 7, corresponding to an average spectral 
luminosity of ~ 5.4 x 10^^ ergs s^^ Hz^ at frequencies from 2.4 x 10^^ to 5.8 x lO^^Hz. The 
X-ray flux from this supernova declined subsequently by ~ 44% in one month (Zimmermann et 
al. 1994; Fransson, Lundqvist & Chevalier 1996). These luminosities are an order of magnitude 
lower than the typical bow-shock emission shown in Figure 3 (except for case 6 which considers 
the ISM, where the two are comparable). X-ray emission has been detected within the first 100 
days of a few other supernovae, such as SN 1980K, SN 19941 and SN 1998bw (see Table 1 in Plan 
1999). SN 1980K had a luminosity of ~ 5 x 10^^ ergs s"^ in the 0.2-4 KeV band on day 44 and 
SN 19941 had a luminosity of 1.6 x 10^* ergs s~^ in the 0.1-2.4 KeV band between day 79 and 
day 85. Both of these luminosities are much fainter than the bow-shock emission we calculated. 
X-rays from the vicinity of SN 1998bw were first detected one day after the explosion, and had a 



- 10- 

luminosity of 5 x 10'*'^ ergs s^^ in the 2-10 KeV band (Pian 1999), corresponding to an average 
spectral luminosity of ~ 2.4 x 10^^ ergs s^^ Hz^^ at frequencies from 4.8 x 10^''' to 2.4 x lO^^Hz. 
This is comparable to the bow-shock emission shown in Figure 3 (except for case 6 with the ISM), 
and so the X-ray emission from SN 1998bw could have resulted from a bow-shock around a jet. 
However, the association of the detected X-ray afterglow with SN1998bw is still not certain due to 
the localization uncertainties (Pian 1999). In both the radio and X-ray regimes, detection of the 
characteristic spectrum and time evolution of the bow shock lightcurve can in principle be used to 
separate it from the emission due to the main supernova shock.Q| 

Finally, we would like to find the maximum distances out to which the bow shock emission 
would be detectable in the radio or X-ray regimes. In the radio regime, we adopt a detection 
threshold of 1 mJy at 10 GHz. In the X-ray regime, we consider the sensitivity of the Chandra 
X-ray Observatory (CXO), corresponding to a flux limit of 2 x 10~^^ ergs cm^^ s^^ in the 0.5-2 
KeV band [= (1.2 - 4.8) x lO^'^Hz] for an integration time of 130 ks (Giacconi et al. 2000). 
Obviously, this fiducial sensitivity can be improved with a longer integration time. The limiting 
distances corresponding to these detection thresholds are listed in Table 1. We find that at 
v = 10^^ Hz the bow shock emission can be detected out to ~ 60 Mpc when the total energy input 
to the jets is equal to 5 x lO^'^ ergs; in this case, the equivalent isotropic energy output of the GRB 
is 10^^ ergs corresponding to case 3 in the calculation. In the X-ray regime, detection is possible 
out to larger distances although there might be a problem of confusion with other sources at the 
arcsecond angular resolution of CXO. 



4. CONCLUSIONS 

Beamed GRBs produce a bow shock in their gaseous environment, which emits non-thermal 
synchrotron radiation extending from radio to X-ray frequencies. We calculated the bow-shock 
luminosity from beamed GRBs in a precursor wind or ISM environments (see Fig. 3), during the 
first few weeks after the explosion. For typical parameters, we find that the bow shock emission 
from a jet with half- angle ~ 5° is visible out to tens of Mpc in the radio and hundreds of Mpc in 
the X-rays (Table 1). We emphasize that the calculated bow shock luminosity is highly sensitive 
to the total hydrodynamic energy carried by the jets and the density of the ambient medium (see 
cases 2, 3 and 6 in Fig. 3). 

The event rate for lines of sight outside the cones of the jet emission is larger by a factor of 
~ 260(Gfe/5°)^^ than for lines of sight which detect the 7-ray burst itself and are constrained to 
pass through the jet. The rate of 7-ray events is estimated to be, ~ 2.5 x 10^^ yr^^ per L^ galaxy 



^Note that in principle, there could also be a contribution from reflected X-rays due to Compton scattering of the 
beamed GRB emission by the ambient gas (Madau, Blandford, & Rees 2000). However, for the explosion parameters 
and photon frequencies we consider here, this component appears to be negligible compared to the bow shock emission. 
Moreover, its decay with time is much faster than that of the bow shock component. 
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for isotropic emission (Wijers et al. 1998). The Virgo cluster has a total luminosity of ~ 430L^ 
(Woods & Loeb 1999) corresponding to a total event rate of ^ 3 x lO~^(06/5°)~^ yr~^ out to a 
distance of ~ 20 Mpc. At larger distances, the mean density of L^ galaxies (Folkes et al. 1999) 
implies that during an observing period Tobs there should be at least one event at a distance of 
D ~ 90 Mpc(robs/10 yr)"^/^(Gb/5°)^/^. These estimates and Table 1 imply that future radio and 
X-ray observations of SN remnants in the local universe can be used to test whether collimated 
GRBs are associated with SNe. 

This work was supported in part by grants from the Israel-US BSF (BSF-9800343) and NSF 
(AST-9900877), and the NASA grant NAG5-7039. 
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Fig. 1. — Geometry of the bow shock. 



case 


time 
(days) 


L^{u = IQio Hz) 
(ergs s-i Hz-i) 


^radio 

(Mpc) 


L^{u = 1.2 X 10^^ Hz) 
(ergs s-i Hz~i) 


(Mpc) 


1 (standard) 


2.7- 13.2 


2.2 X 10^6 - 5.7 X 10^6 


13.6 - 21.8 


2.4 X 10^2 _ 9 4 X io22 


347 - 687 


2 (A = 0.4) 


6.8 - 32.9 


2.9 X 10^6 - 6.5 X 10^6 


15.6 - 23.3 


6.8 X 10^^ - 2.6 X 10^2 


185 - 361 


3 (^52 = 10) 


27.4 - 131.7 


1.1 X 10^7 - 4.2 X lO^^^ 


30.3 - 59.2 


1.3 X 10^2 - 5.2 X 10^2 


255 - 511 


4 (e = 0.5) 


2.7- 13.2 


2.9 X 10^7 - 3.5 X 10^7 


49.2 - 54.1 


1.1 X 10^3 - 4.6 X 10^3 


743 - 1519 


5 (p = 2.2) 


2.7-13.2 


2.7 X 10^5 - 1.4 X 10^6 


4.8 - 10.8 


4.0 X 10^2 - 1.0 X 10^3 


448 - 708 


6 (ISM) 


36.3 - 97.2 


4.2 X 10^5 - 2.6 X 10^6 


5.9 - 14.7 


1.1 X 10^1 - 2.4 X 10^1 


74 - 110 



Table 1: Values of observed time, the luminosity sA, v = 10^*^ Hz, the detection threshold distance 
at V = 10^*^ Hz, the luminosity sA, v = 1.2 x 10^^ Hz, and detection threshold distance at 
z/ = 1.2 X lOi'^ Hz. 
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Fig. 2. — Coordinate system for the calculation of the bow shock emission. 
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Fig. 3. — Lightcurves of the bow shock emission for six different cases (see §3 for details). In each 
case, the start and end times are labeled on the plot. The intermediate times are evenly distributed 
between the start and end times. 



